Radiation reaction and the self-force for a point mass in general relativity 
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A point particle of mass /i moving on a geodesic creates a perturbation /i M , of the spacetime metric 
g°, that diverges at the particle. Simple expressions are given for the singular \ijr part of ft* 1 and its 
quadrupole distortion caused by the spacetime. Subtracting these from /i M leaves a remainder h R 
that is C . The self-force on the particle from its own gravitational field corrects the worldline at 
0(/x) to be a geodesic of g° + h R . For the case that the particle is a small non-rotating black hole, 
an approximate solution to the Einstein equations is given with error of 0(/x 2 ) as /i — > 0. 
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Perturbation analysis provides one approach to under- 
standing the self-force and radiation reaction in general 
relativity. This begins with a background spacetime met- 
ric g° which is a vacuum solution of the Einstein equa- 
tions G ab (g°) — 0. An object of small mass /i then dis- 
turbs the geometry by an amount — O(p) which is 
determined by the perturbed Einstein equations with the 
stress-energy tensor of the object being the source, 

E ab (h») = -%nT ab + 0(p 2 ), (i -» 0. (1) 

The superscript ^ is a reminder that hf 1 is linear in fi. 
The linear differential operator E ab is defined by 

E ab (h) = -^hd, (2) 
og c d 

which evaluates to (!]] 

2E ab (h) = V 2 h ab + V a V b h - 2V (o y c h b)c 

+2R a c b d h cd + g ab {V c V d h cd -V 2 h). (3) 

With a solution of Eq. (|l|) it follows that 

G ab (g° + hn =8ttT o6 + cV). (4) 

An intcgrability condition for Eq. (|l|) results from the 
Bianchi identity for g° + and requires that T be con- 
served in the background geometry up to 0(^ 2 ). 

Formally, perturbation analysis at the second order is 
no more difficult than at the first. But the integrability 
condition for the second order equations is that T be con- 
served not in the background geometry, but in the first 
order perturbed geometry. Thus, before solving the sec- 
ond order equations, it is necessary to change the stress- 
energy tensor in a way which is dependent upon the first 
order metric perturbations. This correction to T is said 
to result from the "self-force" on the particle from its own 
gravitational field and includes the dissipative effects of 
what is often referred to as "radiation reaction" as well 
as other nonlinear aspects of general relativity. 

To focus on those details of the self-force which are in- 
dependent of the object's structure we restrict the object 
to be a point particle with no spin angular momentum 



or other internal structure. The integrability condition 
at the first order then implies that the worldline V of the 
particle is nearly a geodesic of g , with an acceleration 
of only 0(/i). But, the integrability condition at the sec- 
ond order presents a difficulty. The particle is to move 
along a geodesic of g° + h^, but scales as /x/r near the 
particle and is not differentiable on T. 

Mino, et al. j| and Quinn and Wald |J resolve this 
difficulty with a Green's function approach to Eq. (Q). 
The formal Hadamard expansion of the Green's function 
near the worldline of the particle identifies the "instan- 
taneous" and "tail" parts of And, Mino et al. use 
a matched asymptotic expansion to show that the parti- 
cle moves along a geodesic of g° + h% ail ; by construction 
^taii ^ s differentiable on the worldline as required by the 
geodesic equation. However, their analysis provides no 
simple method for determining this tail part. 

Alternatively, we resolve the difficulty by finding the 
source part of the metric perturbation h s , which con- 
sists of the singular ji/r part plus its quadrupole distor- 
tion caused by the background geometry. Eqs. (|lO|)-(|T^) 
give a simple expression for h s . Then we show that the 
remainder, h R = — h s , is C 1 and, using matched 
asymptotic expansions, that the 0(/i) effect of the self- 
force adjusts the worldline of the particle to be a geodesic 
r' of g° + h R . The consistency of our matched asymptotic 
expansions with those of Ref. || imply that h R must be 
equivalent to the "tail" part of the metric perturbation 
from the Green's function, up to a gauge transformation 
and terms of 0(fxr 2 ), which do not effect the 0(/x) cor- 
rection to the worldline. 

The source field h s is best described with coordinates 
in which the background geometry looks as flat as pos- 
sible near the geodesic T. A normal coordinate system, 
x a = (t,x,y,z), can be found [jj] where, on T, the met- 
ric and its first derivatives match the Minkowski metric, 
and the coordinate t measures the proper time. Normal 
coordinates for a geodesic are not unique, and we use 
particular coordinates introduced by Thorne and Hartle 
H in their discussion of external multipole moments of 
a vacuum solution of the Einstein equations where 

gib = Vab + iH ab + 0(r 3 /^ 3 ), r/K 0, (5) 
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with 

2H ab dx a dx b — ~£ijX l x^ (dt 2 + Skidx k dx l ) 

+ je kpq B q l x p x t dtdx k . (6) 

And, £ and B are spatial, symmetric, tracefree and re- 
lated to the Riemann tensor evaluated on T by £ij = Rutj 
and Bij = ti pq R pq jt / '2; and, TZ is a representative length 
scale of the background geometry — the smallest of the ra- 
dius of curvature, the scale of inhomogeneities, and the 
time scale for changes in curvature along T, then £y and 
B q i are 0(1/1Z 2 ); also (r,Q,<f>) are defined in the usual 
way in terms of (x, y, z); the indices i, j, k, . . . are spatial 
and raised and lowered with Sij. 

If a small non-rotating black hole moves along T, then 
its geometry is perturbed by tidal forces, 

„pcrt _ Schw , ; f 7 \ 
9 a b - 9ab +2llab {() 

through terms of 0(r 2 /TZ 2 ), where ih is a solution of 

^(2*) - (8) 

with the boundary conditions that the perturbation be 
well behaved on the event horizon and that ih — > 2H in 
the buffer region Q], where /1 <C r <C TZ. Both 2H and 2/1 
consist of £ = 2 tensor harmonics in the Regge- Wheeler 
gauge ||; the angular dependence is through x l x^ £ij and 
tkpqB q iX p x l . For r < 2/1 is governed by a wave equa- 
tion with a potential barrier. In the time independent 
limit this admits an analytic solution Q 

2 h ab dx a dx b = -£ ij x i x j [(1 - 2[i/r) 2 dt 2 + dr 2 
+ (r 2 - 2n 2 )(d0 2 + sin 2 6d<j, 2 )} 

+ ^ e^BO^x 1 (1 - 2/i/r)dt (9) 

which is well behaved on the event horizon and matches 
2 if when fi <C r. See Ref. [|| for the case of a small black 
hole in the vicinity of a larger black hole. 

Time dependence of 2H induces a quadrupole mo- 
ment on the black hole, but the resulting acceleration 
of the world line is smaller than 0(fj,/TZ 2 ). Generally, 
the timescale of 2H is TZ and corresponds to a low fre- 
quency for the black hole, cjfi = 0(fi/TZ) <C 1. And two 
independent solutions for the metric perturbation 2/1 are 
standing waves very near r = 2/i but behave as r 2 and 
1/r 3 , for /i « r « K. The proper solution is a travel- 
ing wave into the hole created by a linear combination of 
these two independent solutions having comparable mag- 
nitudes when r ~ 2\i. Thus, for r>/i this linear combi- 
nation is approximately given by Eq. (Q) , which scales as 
r 2 , plus a 1/r 3 contribution from the induced quadrupole 
moment, T ab = 0(p 5 /TZ 2 ), stemming from the time de- 
pendence. This contribution to the quadrupole field cou- 
ples to the background octupole field and accelerates the 
worldline by ~ £ i ab l ab /fi = 0(fi 4 /TZ 5 ) which is too 
small to be important in this analysis. 



In the buffer region, where fi <C r <C TZ, the geometry 
of a point particle moving through the background should 
be equally well described either by the background met- 
ric perturbed by /1, or by the leading ji/r terms of the 
Schwarzschild metric perturbed by weak tidal forces. In 
this region, then, the background metric perturbation 
is approximately the part of g pert which is linear in /i; and 
this part is the source field, 

h s = ^ + 2 h^, (10) 

where 

h^ b dx a dx b = 2fi/r(dt 2 + dr 2 ) (11) 
is the fi/r part of the Schwarzschild metric, and 

2 h^ b dx a dx b = ^£ tjX l x 3 dt 2 - ^-e kpq B\x p x l dtdx k 

(12) 

is theur/TZ 2 part of 2^ from Eq. (^|). A split of from 
Eq. ([j]), into two pieces, 

h» = h s + h R , (13) 

reveals just how accurately h approximates /i M by con- 
sideration of the remainder h R . From Eq. (|l]) 

E ab {h R ) = ~E ab {h s ) - 8nT ab . (14) 

And, direct evaluation shows that 

E ab {h s ) + 8nT ab = 0(n/TZ 3 ), r -» 0, (15) 

for a point particle stress tensor, and the source of Eq. 
( |T4| ) is finite but not continuous at r = 0. This last result 
is understandable — in an expansion of E ab (h ) in powers 
of fi and 1/TZ, all of the [i/TZ 2 terms would also appear in 
a similar expansion of Ef b (2J1). And this latter expan- 
sion is zero from the definition of 2 h. Thus, the source of 
Eq. d||) is 0{^/TZ 3 ) as r -> 0. 

The solution of Eq. ( pT| ) for h R with reasonable bound- 
ary conditions is C . If the metric is sufficiently differen- 
tiable and there are no unreasonable boundary conditions 
at large distances, then h R is differentiable away from the 
geodesic T. For if it were not, then the discontinuities in 
the derivatives would propagate along the characteristics 
of the hyperbolic operator E ab and would have originated 
either on T or on some boundary; we consider such dis- 
continuities emanating from a boundary to be unreason- 
able boundary conditions. And, in a neighborhood of 
r the geometry can be smoothly mapped to flat space- 
time with the operator E ab (in the Lorentz gauge) being 
smoothly mapped to the flat spacetime wave operator 
which, when integrated twice, smoothes a slowly chang- 
ing and finite but discontinuous (on V) source into a C 1 
solution. Thus the difference between N 1 and h s is a C 1 
tensor field h R . 
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The split of W 1 into h s and h R contains some arbitrari- 
ness. Any piece of 0(^,r 2 /TZ 3 ) or with higher powers of 
r/TZ can be moved between h and h R without affecting 
either Eq. ([il]), the differentiability of h R , or the O(fi) 
effect of the self- force which changes the worldline to be a 
geodesic V of g° + h R . Furthermore, a gauge transforma- 
tion y a = x a + t; a for any £ a that is 0(r 3 /TZ 2 ) gives a new 
normal coordinate system; and we state without details 
that the corresponding change in h s is only 0(fir 2 /TZ 3 ) 
with vanishing derivatives on the worldline. Thus T' is 
independent of the normal coordinate system in use. 

But, it is not sufficient to have h s just consist of oh' 1 . 
If 2^i M were not included, then E a t,(h s ) + %irT a i, would be 
singular ~ fi/rTZ 2 as r — > 0. The resulting h R would not 
be differentiable on the worldline, and some version of 
averaging around the particle would be required to make 
sense of the effects of the self-force. Thus, it is necessary 
to include 2^ as part of h s . 

Above, we mentioned that the worldline of a small par- 
ticle through the background is a geodesic of g°+h R when 
the O(n) corrections are included. We now justify this 
statement by replacing the particle with a small, non- 
rotating black hole and considering a sequence of metrics 
g(fj-) which are solutions of the vacuum Einstein equa- 
tions, and parameterized by \i with g(0) = g°. Our focus 
is on the 0{fi) behaviors of g(/i) and the worldline in the 
limits /it — ► and the worldline approaching T. Match- 
ing asymptotic expansions in the buffer region provides 
both the correction to the worldline as well as an approx- 
imation to gin), Eq. ( |l9| ) below, which is uniformly valid 
with an error of 0(p? /TZ 2 ), as fj, — > 0. 

This matching relies upon a choice of coordinates sim- 
ilar to that of Eq. (||) but for the metric g° + h R being 
expanded about its geodesic T'; these "primed" normal 
coordinates differ from the original by 0(/i). In this dis- 
cussion of matching, r is now associated with V and the 
primed coordinates. 

In the buffer region g(fj,) is best described in a fashion 
introduced by Thorne and Hartle Q as a sum of terms 
of positive powers of the small numbers \i/r and r/TZ, 

g(n) ~ rf & & % ff & ••• 

& fi/r & n/TZ & \ir/TZ 2 & ••• , , 
& fi 2 /r 2 & fi 2 /rTZ & fi 2 /7Z 2 & •■• [W) 
& 

where & means "and a term of the form . . ." In this 
tableau, part of any term can be moved diagonally, up 
and to the left, to be absorbed into a dominating term 
with the same r behavior. 

Just outside the buffer region, where fi <C r < 7Z, g(n) 
is approximately the background geometry perturbed by 
fi. The top row of the tableau consists of the expansion 
of g° about V in powers of r/TZ along with some parts 
of a similar expansion of h R . The fi/TZ and fir/TZ 2 terms 
from h R are absorbed into the first two elements of the 
top row, which are then displayed as rf and zero because 
the primed coordinates are normal. The r 2 /TZ 2 term of 



the tableau is set to 2H', a tensor whose components in 
the primed coordinates are the same as those of 2H in the 
original coordinates. The difference between 2H' and 2H 
is of 0([ir 2 /TZ 3 ) and is subtracted from the fourth term 
in the second row of the tableau. The rest of the first 
row consists precisely of the corresponding terms in the 
expansion of g° . And the remainder of the expansion of 
h R , from the [ir 2 /TZ 3 term outward, contribute to the 
corresponding terms of the second row. Thus the entire 
first row along with the fourth and greater terms of the 
second row sums to g° + h . 

Just inside the buffer region, where 2^ < r < R, 
is approximately the Schwarzschild geometry per- 
turbed by background tidal forces. The first column of 
the tableau, containing no TZ, is an expansion of the 
Schwarzschild geometry in powers of /i/r. The second 
column, linear in 1 /TZ, sums to a dipole perturbation of 
the Schwarzschild geometry. But, the top element of the 
second column is zero, so all elements of the second col- 
umn are zero. The top term in the third column, 2H', 
when added to the rest of the third column gives 
the quadrupole perturbation of the black hole caused by 
tidal forces. Thus, the first three elements of the top row 
determine the entire first three columns of this tableau 
by the expansions of g Schw and 2^ in powers of (j,/r. 

Now, g° + h is an accurate approximation of g(fj) 
when \x r, and g Schw -f 2 h' is an accurate approxima- 
tion when r -C TZ and with g Schw centered on V. These 
approximations overlap in the buffer region where 

9ab + hab - Vab + 2Kb + 0{r 3 /TZ 3 ), (17) 

and 

9 S ab hw + 2Kb = Vab + 2Kb + 0(/i/r) (18) 

match asymptotically. In the restricted region /i/r 
r 2 /TZ 2 < 1, the displayed term 2 H' = 0(r 2 /TZ 2 ) is 
small yet much larger than either of the remainder terms, 
0{r 3 /TZ 3 ) and 0(n/r), as fi/TZ -> 0. This is the hallmark 
of matched asymptotic expansions. If the worldline of 
the black hole were not a geodesic of g + h R then Eq. 
( |l7| ) would necessarily contain a term of 0(r/TZ). But 
an explicit 0(r/TZ) term in Eq. ( |l8| ) would be the dom- 
inant term of a dipole perturbation, and such a dipole 
vacuum perturbation of the Schwarzschild geometry is al- 
ways removable by a gauge transformation g. Thus, this 
asymptotic matching is only successful when the world- 
line r' is a geodesic of g° + h R up to an acceleration of 
0(fi 2 /TZ 3 ) in the limit that fi/TZ -» 0. 

A concise description of this matched geometry is 

9ab{^) = {gib + h R b ) + (g s a t lw + 2h' ab ) 

- Kb + 2Kb) + O^/TZ 2 ), n/TZ - 0. (19) 

For r <^ TZ, the first and third terms on the right nearly 
cancel and give <?(/i) « g Schw +2h' , the first three columns 
of the tableau. For fi r the second and third terms on 
the right yield h s ' + 0(/j, 2 /r 2 ). And g(fi) « g° + h R + h s ', 
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the top two rows of the tableau. The error in this approx- 
imation comes from the n 2 r/7Z 3 term in the tableau, the 
dominant term not fixed in the above discussion. This 
term and the rest of the third row sum to a contribution 
of 0(fi 2 /1Z 2 ). Thus, the matched asymptotic expansions 
give an approximation to in Eq. (1^) which is uni- 
formly valid with error of 0(fi 2 /7Z 2 ) as \i — > 0. 

An application of this approach, in conjunction with 
Fourier-harmonic decomposition, determines the 0(/x) 
corrections to geodesic motion for a small non-rotating 
mass in orbit about a much larger non-rotating black 
hole. First, Eq. (|l|) is solved for using the usual met- 
ric perturbation analysis of the Regge- Wheeler formal- 
ism ||. This involves decomposing T into its Fourier- 
harmonic modes; then, the inhomogeneous Regge- 
Wheeler or Zerilli equation Q is integrated numerically 
to determine the radial dependence of the modes of h^. 
And, given the appropriate coordinate transformation, 
the components of h s can be transformed from Eqs. 
( pd| ) and ([l2|) to the usual Schwarzschild coordinates 
and then numerically decomposed into their Fourier- 
harmonic modes. Now, h R and its derivatives can be 
constructed as the sum over modes of the difference be- 
tween W and h s . And, the 0(/i) effect on the worldline 
of the small mass can be calculated as a change from a 
geodesic of g° to a geodesic of g° + h R . 

This application depends upon the transforma- 
tion between the normal coordinates of Eq. (H) and 
Schwarzschild coordinates. Manasse and Misner B give 
a prescription for finding Fermi-normal coordinates for 
any geodesic; and Zhang || gives a gauge transformation 
from these to the coordinates of Eq. (H) . But Eq. (||) only 
determines the normal coordinates near the worldline and 
up to terms of 0(r 4 /7Z 3 ). So, the normal coordinates 
must be extended to cover the Schwarzschild manifold 
for all radii in the vicinity of the orbit to make the mode 
decomposition of h s possible. Fortunately, the details of 
this extension are not important because h R is known 
to be differentiable and, therefore, easily describable in 
terms of a sum over its modes. While the amplitude of 
each individual mode of h does depend upon the ex- 
tension of the coordinates away from the worldline, the 
reconstruction of h R near the worldline by the sum over 
modes is independent of this extension. 

This analysis shows that simply removing the diver- 
gent monopole part from the metric perturbation 
/i M leaves a nondifferentiable remainder. But, if the 
quadrupole distortion of the monopole is also removed 
from h^, then the remainder, h R — U L — W l — 2 h^, is 



differentiable and suitable for calculating 0(/i) effects on 
the worldline. 

The octupole term 3/1^ can also be removed from . 
Thorne and Hartle 0] extend the coordinates in Eq. (|5|) 
to include 0{r 3 /1Z 3 ) terms explicitly, their Eqs. (Al) and 
(A2). The time independent solution of E^ hw ( 3 h) = 
which is well behaved on the event horizon and properly 
matches the 0(r 3 /1Z 3 ) terms of Ref. [Q is 



3 h ab dx a dx b = -1^iVi*[(1 - ^) 2 (1 

4^ 3 



/'■■ 



dt 2 



(1 _ £) df 2 + ( r 2 _ 2flr + ^L)(^2 + sin 2 
T 07* 



<h-^)a-^E) dtdx \ (20) 

r or 



And, 3-ff and 3/1^ are the r 3 /1Z 3 and /ir 2 /1Z 3 parts of 
3/1. If 3/1^ is also removed from h 1 * then the remain- 
der is C 2 , the matching is extended through the 1/1Z 3 
terms, but the overall error of the matched geometry is 
still 0(p 2 /7Z 2 ) and there is no effect on the worldline. 

Interesting discussions with Bernard Whiting, Eric 
Poisson, and Yasushi Mino greatly clarified many aspects 
of this letter. This work was supported in part by NASA 
grant NAGW-4864 with the University of Florida. 
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